COLO-HEP-465 
May 2001. 

Four Flavor Finite Temperature Phase Transition with 
HYP Action: Where is the First Order Phase Transition 

Line? 

§ ! Anna Hasenfratz^ and Francesco Knechtli* 

<n : 

Physics Department, University of Colorado, 

2 . 

•"■> ' Boulder, CO 80309 USA 

o : 

Abstract 

<N ■ 

| We study the finite temperature phase transition of four flavor staggered fermions with 

(N ■ 

hypercubic fat link actions on N t = 4 and N t = 6 temporal lattices. Our fat links are 

in : 

constructed with hypercubic blocking (HYP) and therefore are very compact. We present 
a new algorithm for simulating fermions coupled to HYP fat links. The algorithm has a 
simple form based on the standard overrelaxation and heatbath updatings for the pure 
Ph! gauge action. We observe that as we increase the smoothness of the gauge fields by 

changing the parameters of the blocking the very pronounced first order phase transition of 
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the thin link action becomes weaker and moves to physically uninteresting values of the 



gauge coupling. With our smoothest HYP action we do not find any indication of a phase 

transition in the accessible temperature range on the N t — 4 or 6 lattices even at quark 
masses close to the physical light quark mass. We argue that the observed difference in the 
phase diagram is due to the improved flavor symmetry of the smeared link actions. 
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I. INTRODUCTION 



Understanding the finite temperature phase structure of QCD has always been at the 
center of lattice studies. Recently all efforts have been focused on two and 2+1 flavor 
investigations 0, hardly any work looked at four flavor QCD |2|, |3|, £|, |5], || 0, 0. This is 
not surprising, since early numerical results confirmed the theoretical expectations of the 
four flavor QCD phase diagram and computer resources were quickly directed towards the 
physically more interesting two and 2+1 flavor cases. 

Based on universality Pisarski and Wilczek argued that QCD with rif > 3 massless 
quark flavors has a first order chiral phase transition ||5J. At finite quark mass the n\ — 1 
degenerate massive Goldstone bosons can destroy the phase transition. Above some critical 
mass value m max one expects only a crossover. Numerical results using thin link staggered 
fermions confirmed the first order phase transition and determined the critical endpoint on 
N t = 4, 6 and 8 lattices. Since the pure gauge SU(3) theory has a first order deconfining 
phase transition, at very large quark masses one might again observe a first order transition 
but numerical simulations have not explored that region. It is only a small shadow on the 
otherwise satisfying picture that the numerically observed critical mass, m max , does not 
scale as the temporal lattice size increases from four to six. This small "imperfection" can 
be explained by scaling violations at large lattice spacing 0. 

Just because numerical simulations seemingly agree with theoretical expectations does 
not mean that we have the correct physical picture. The Pisarski- Wilczek argument is 
based on the SU{rif) x SU{rif) symmetry of the fermion action and the existence of nj — 1 
degenerate Goldstone bosons in the chirally broken phase. Staggered fermions, on the other 
hand, break flavor symmetry. There is only a remnant U(l) symmetry and a single true 
Goldstone boson in staggered simulations. The other, "would-be" Goldstone bosons are 
massive pseudoscalar particles even at vanishing quark mass. The full chiral symmetry is 
recovered only in the continuum limit. At the lattice spacing corresponding to N t = 4 critical 
temperature the flavor symmetry violation of the thin link staggered action is substantial, 
the Pisarski- Wilczek scenario is not applicable. 

Chiral symmetry breaking in QCD is closely related to the topological structure of the 
vacuum. Recently it was shown that the topological susceptibility measured on two- or 
four-flavor thin link staggered fermion configurations at lattice spacing a ~ 0.17fm does not 
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follow the theoretically predicted chiral behavior at small quark masses [pTOfl . Rather, the 
topological susceptibility appears to be consistent with the quenched value independently 
of the quark mass, implying that the thin link staggered fermions at that lattice spacing do 
not act like two or four degenerate flavors. The corresponding vacuum, at least regarding 
instantons, is closer to the quenched than to the dynamical vacuum. One has to reduce 
the lattice spacing to a ~ 0.1 fm to get acceptable, though still slightly high, values for the 
topological susceptibility. If instantons are responsible for chiral symmetry breaking, the 
chiral finite temperature phase transition should show similar lattice artifacts. The critical 
temperature of the chiral restoring phase transition is expected to be T c ~ 150 — 170 MeV 
||. A lattice spacing a ~ 0.17 fm then corresponds to a critical system with temporal lattice 
extension N t = l/(Ta) ~ 7 — 8, a lattice spacing a ~ 0.1 fm corresponds to a system with 
N t ~ 12 — 14. Therefore, if the indications from topology apply to the chiral phase transition, 
one has to use very large lattices, N t > 8, to study the finite temperature phase transition 
with thin link staggered fermions. 

There is no strong dependence on the number of fermion flavors in the above argument, 
assuming that the critical temperature does not depend strongly on n/. Instanton model 
calculations indicate that this might not be the case. Results obtained using the Interacting 
Instanton Liquid Model (IILM) suggests a strong flavor dependence |ITJ, |12] . For rif = 2 
flavors the model predicts a second order transition at T c « 150MeV in the chiral limit and 
only crossover at finite quark mass. With rif = 3 flavors the chiral transition becomes first 
order and it occurs at considerably lower temperatures, T c rs lOOMeV. With n/ = 5 flavors 
the IILM predicts that there is no spontaneously broken phase in the massless chiral limit, 
that phase exists only at non- vanishing quark mass values. The rif = 4 case appears to be 
borderline. If a non-zero chiral condensate exists, it is very small and occurs only at very low 
temperatures. These considerations indicate that in order to reproduce the physical phase 
diagram on the lattice the instanton content of the vacuum should be correctly reproduced 
by the lattice action. 

Flavor symmetry can be considerably improved by coupling the fermions to fat links 



14, 15fl. Improved flavor symmetry leads to improved topology. The first results of Ref. 



T0|| indicate that fat link actions can reproduce the quark mass dependence of the topological 



susceptibility even at lattice spacing a ~ 0.17fm. It is possible that even coarser lattices 
can be used with fat link fermions though it is unlikely that a much larger lattice spacing is 
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allowed because instantons fall through the lattice when their radius p is comparable to the 
lattice spacing, p/a ~ 1. Since the average instanton radius is p ~ 0.3 fm, configurations 
with a > 0.3 fm will not have correct chiral behavior. Because of that dynamical lattices 
with temporal extension N t — 4 are probably too coarse to reproduce the correct continuum 
chiral behavior even with chiral fermions, N t = 6 — 8 might be acceptable with chiral or 
almost chiral fermions. 

In a recent publication we proposed an algorithm to simulate dynamical fat link actions 
TEfl . In this paper we use a modified version of that method to study the phase diagram of 



four flavor staggered fermions on Nt — 4 and N t = 6 temporal lattices. We use a new type 
of fat link action created with a hypercubic block (HYP) transformation [16|]. The HYP 
blocking mixes gauge links within hypercubes attached to the original link only so that 
lattice artifacts due to extended smearing are minimized. The HYP fat links are so compact 
that the static potential is indistinguishable from the thin link one at distances r/a > 2, 
yet flavor symmetry with staggered quarks is improved by about an order of magnitude. 
Our results for the finite temperature phase diagram with HYP staggered fermions deviate 
significantly from the thin link action predictions. At N t = 4 and constant physical quark 
mass values the strongly first order phase transition observed with thin links weakens and 
moves deep into the strong coupling region as fattening is introduced. It appears that the 
N t — 4 first order phase transition with thin link staggered fermions is more a lattice artifact 
than a physical phase transition. With our smoothest HYP action we do not see any sign of 
a first order phase transition even at quark masses comparable to the physical light quark 
mass, neither on N t — 4 nor on N t — 6 lattices. Because of the large lattice spacing at 
N t — 4 even a chirally symmetric action might not be reliable but the N t = 6 results should 
be closer to the continuum behavior. We believe that the observed difference between the 
HYP and thin link actions is due to the 15 near-degenerate Goldstone bosons of the fat link 
HYP action which can change the phase diagram. 

There are several groups doing calculations using dynamical fat link actions though with 
fat links that are less smooth than the HYP fat links |17|, [L8|. Unfortunately they did 
not publish results about the finite temperature phase transition with four flavors or on 
the topological susceptibility with either two or four flavors using fat link actions. We feel 
that the topological susceptibility is one of the best indicators to distinguish quenched and 
dynamical configurations and verify the effect of the sea quarks on the vacuum. Since Ref. 
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T0|| has only preliminary results using a different fat link action than considered here, we are 



presently investigating the topological susceptibility of the HYP action at different lattice 
spacings and quark masses. Early results show consistency with the theoretically predicted 
chiral behavior at a lattice spacing a ~ 0.17fm and at a quark mass corresponding to 
a Goldstone pion mass m^ro — 2.0 . Detailed results will be presented in a forthcoming 
publication 

This article is organized as follows. In Sect. 2 we present our new algorithm to simulate 
fat link actions. It is based on the standard overrelaxation and heatbath updatings for the 
pure gauge action. We describe the actions used in this study and the performance of the 
algorithm. In Sect. 3 we present our finite temperature results for four flavors of staggered 
fermions. We show how the phase diagram at N t = 4 is changed by smoothing the gauge 
fields. With our smoothest HYP action we study the phase diagram when the quark mass 
is lowered and the spatial volume is increased. We also simulate the HYP action on finer 
lattices at N t — 6. There is no sign of a first order phase transition even with very light 
quarks in the physically relevant coupling range. In Sect. 4 we conclude with a summary. 



II. SIMULATIONS 



A. Algorithm 



Dynamical simulations of fat link actions are difficult because the fermions couple to a 
complicated, extended, and, in case of projected fat links, non-linear combination of the 



gauge links. In |15[ we considered fat links created by several levels of projected APE 
blocking. We introduced an auxiliary gauge field for each blocking level and coupled the 
fermions to the last level of auxiliary gauge fields. 

In this article we present a new action that has no auxiliary gauge fields and can be 
simulated in a simpler way. Our new fat link action is 

S = -^J2 ReTr ( U p)- tr HM*(V)M(V)}, (1) 
3 p 

where U p is the plaquette product of the thin gauge links £Zj i(U , M is the fermionic matrix and 
V i)fJi is the fat link that couples to the fermions. We denote by Tr the trace over SU(3) color 
whereas tr means the trace over space-time indices i, directions fi, spin and color. If the fat 
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links are constructed in a deterministic way from the thin gauge fields, they are no longer 
dynamical variables and an ergodic update of the U fields will correctly simulate the system. 
Fat links can even be constructed by iterating the block transformation several times and 
projecting back onto SU(3) after each blocking step. The general algorithm to simulate the 
action of eq(|I]) is based on a two-step decomposition of the action. First, a subset of the 
thin links are updated with the standard microcanonical overrelaxation [^, ^IJ or Cabibbo- 
Marinari heatbath [^] algorithms for the pure gauge action Sg(U) = —^J2 p R e Tr(U p ). 
Next, the deterministically constructed fat links are updated and this change is accepted 
with probability 

Pacc(V',V) = mm{l,exp[-S F (V) + S F (V)]} , (2) 

where V denotes the new fat link configuration and S F (V) = -tr \n[M\V)M{V)} is the 
fermionic action. The sequence of thin link updates must be carefully chosen to satisfy 
detailed balance with respect to the pure gauge action Sq- The proof that this algorithm 
satisfies detailed balance with respect to the full action S = So + S F follows closely the 
proof given in ||23|| for a similar two-step decomposition of the action. 

In the practical implementation of the algorithm we choose the subset of thin links to be 
updated differently for the overrelaxation and for the heatbath steps. For the overrelaxation 
we reflect all the links within some finite block of the lattice. The location of the block 
is chosen randomly but its dimensions are fixed. We choose with probability 1/2 a given 
sequence of reflections within this block or with equal probability the reversed sequence. The 
sequence has to be reversed with respect to the direction and location of the thin links and 
with respect to the index of the SU(2) subgroup used in the reflection step. For the heatbath 
we choose the subset of thin links randomly, i.e. we choose a random direction, a random 
parity, a set of random sites and a random sequence of SU(2) subgroups. The probability 
of generating a given sequence of thin link updates or the reversed sequence is then the 
same, both for overrelaxation and heatbath. From this it follows that detailed balance with 
respect to the pure gauge action Sq is satisfied for both updatings. We chose to update 
a contiguous block of links with overrelaxation in order to propagate a change through (a 
portion of) the lattice. We observed that a random sequence of overrelaxed reflections is 
less efficient, especially for changing the topological charge of the configurations. This can 
be understood by thinking of instantons as extended objects. One has to change an entire 
region of the lattice to destroy or to create them. 
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The acceptance probability eq(§) contains the ratio of fermionic determinants 
&et[M\V')M(V')}/ &et[M\V)M(V)} which we evaluate by the same method as described 
in |15|. To make this article self-contained we summarize the steps in the following. First, 
we remove the most ultraviolet part of the fermion matrix by decomposing it as 

M(V) = M r (V)A(V) with (3) 

A(V) = exp [a 4 D\V) + a 2 D 2 {V)] , (4) 

where D is the kinetic part of the fermion matrix and the parameters a± and a 2 are arbitrary 
but real. This way we achieve that an effective gauge action 

S eff (V) = -2a A Retr[D 4 (V)}-2a 2 Retr[D 2 (V)} (5) 

is removed from the determinant. The acceptance probability eq(|2|) can be approximated 
by the stochastic estimator 

P' acc {VW) = min{l,exp{5' e// (y) - S eff (V) + £ f [M}(V')M r (V) - M}(V)M r (V)} £}} , 

(6) 

where the vector £ is generated according to the probability distribution 

P(0 oc exp{-eM},(V')M r (V')t} . (7) 
The parameters a 2 and «4 can be optimized to maximize the acceptance rate. We keep the 



same choice a 2 = —0.18 and a A = —0.006 as in 15 



B. Actions 



In this study of the QCD thermodynamics we use four flavors of staggered fermions and 
compare results obtained with three different fermionic actions. Besides the standard action 
with thin link gauge connections, we simulate two fat link actions of the form eq(|l|): the 
APE1 action and the HYP1 action. The simple form of the action and of the algorithm 
based on overrelaxation and heatbath updatings with respect to the gauge action allows us 
to simulate fat link action with any fattening procedure. No further improvement of the 
gauge and fermionic actions is used here. 
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Action/^ (it) 


Thin 


APE1 


HYP1 




0.594(25) 


0.212(22) 


0.086(14) 




0.72(6) 


0.35(5) 


0.150(24) 



Table I: The parameter 82 measuring flavor symmetry violations as computed in the quenched 
approximation at j3 = 5.7 for the different quark actions that we use in this study. 

In the APE1 action the fat links are constructed using only one level of APE smearing 
with APE parameter a = 0.7 ||24||. The fat link is defined as 



QJ ^ 

Vi !tl = Proj S u(3)[(l-a)U i>fi + - ^ U itU U i+u ^ i+ ~ lvi , 



,.U, 



(8) 



where the projection onto SU(3) is deterministic. In the HYP1 action the fat links are 



constructed using one level of hypercubic (HYP) blocking and we refer to Ref. |16| for the 
precise definition and the parameter choice of the HYP fattening. The HYP fat links mix 
thin links within the hypercubes attached to the original link only. The artifacts due to 
extended smearing, e.g. in the static potential at short distances, are minimized with the 
HYP blocking. At the same time the improvement of chiral or flavor symmetry is substantial. 



In |16|] it has also been demonstrated that one level of HYP blocking remains effective as the 
lattice spacing is decreased. Finally, the APE smearing can be considered as a particular 
choice of the parameters of the HYP smearing. 

To compare the flavor symmetry violations of the standard, APE1 and HYP1 actions we 
calculated the quenched pion spectrum on an ensemble of j3 = 5.7, 8 3 x 24 pure gauge con- 



figurations with lattice spacing a ~ 0.17 fm. As we demonstrated in fll5[ the quenched result 



is a good indication of the dynamical flavor symmetry violations. We use the parameter [13 



do 



ml — m, 



a 



m. 



(9) 

»p m G 

at m n /m p = 0.55 to quantify the difference between the mass of the true Goldstone pion 
me and the mass of the pseudoscalar particles (or pions) m n . The results are summarized 
in table |. We compute the parameter 8 2 for the two lightest non-Goldstone pions, 71^5 
with flavor structure 7^75, and n it j with flavor structure 7*7^. For a flavor symmetric action 
^2 = 0. The flavor symmetry violation of the thin link action is reduced by a factor 3 with 
the APE1 action and by a factor 6 with the HYP1 action at this lattice spacing. 
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Figure 1: The mass renormalization in the quenched approximation at (3 = 5.7 for different quark 



actions. The APE3 action is the fat link action of Ref. |15|. The three lowest quark mass values 
with thin link action are from Ref. VQl 



In figure ffl we plot the square of the Goldstone pion mass mg in units of the Sommer 



scale r |26[ as function of the bare quark mass am q . In addition to the three valence actions 
used in this study we plot the data for the APE3 action (3 levels of APE smearing with 
a = 0.7 and projection parameter A = 500) of Ref. [|15|]. The data for the three lowest quark 
masses with thin link action come from Ref. E3| and are obtained on (3 = 5.7, 16 3 x 32 



lattices. The mass renormalization after one level of HYP blocking is the same as after 3 
levels of APE blocking. We will use these quenched results to approximately match the 
physical quark masses between dynamical simulations of four flavors of HYP and thin link 
staggered fermions. 
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C. Performance 



We would like to compare the performance of the algorithm described in section II A 
for the simulation of the HYP1 action with the standard HMC algorithm for the thin link 
action. We would like to emphasize that the implementation of the algorithm for the HYP1 
action is very simple and can be easily parallelized. We also optimized our parallelized code 
reducing considerably the simulation time cost. 

On 8 3 x 24 lattices we simulated the HYP1 action with gauge coupling (3 = 5.2 and 
quark mass am q = 0.1. An approximate matching of the lattice spacing and the physical 
quark mass can be achieved by simulating the thin link action with parameters f3 = 5.2 and 
am q = 0.06. In an attempt to make a reliable comparison of the time costs, we estimated the 
autocorrelation times for simple observables like the plaquette and the chiral condensate < 
iptp >. We found that measurements with thin link action separated by one HMC trajectory 
of unit time length have comparable autocorrelations as measurements with HYP1 action 
separated by 160 overrelaxation steps and 80 heatbath steps. Each overrelaxation step 
updates a block of 128 links, each heatbath step updates 200 links with an acceptance of 
about 20% for both updatings. In these physical time units the simulation of the HYP1 
action is a factor of 7 more expensive than the thin link action. 

Interestingly neither the heatbath nor the overrelaxation algorithm looses efficiency as 
the quark mass is lowered. With a combination of the two updates one can easily simulate 
systems at very small quark masses that would be impractical if not impossible to simulate 
with standard fermionic algorithms. We simulated on 8 3 x 24 lattices the HYP1 action 
with parameters (5 = 5.2, am, = 0.04 and the thin link action with parameters (3 = 5.2, 
am q = 0.024, with approximately matched physical quark masses. The time costs are now 
only a factor 3 higher for the HYP1 action as compared to the thin link action, due to 
the considerable increase of conjugate gradient steps needed in the inversion of the fermion 
matrix with thin links. 

At a fixed lattice spacing the necessary number of overrelaxation and heatbath updating 
steps between independent configurations scales with the volume. Since each updating 
requires the evaluation of the fermionic action, this gives a volume square dependence for 
the cost of the algorithm. On the other hand the number of links that can be updated 
with the overrelaxation or heatbath algorithms scales with the lattice spacing, the necessary 
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number of updating steps for a fixed physical volume is independent of the lattice spacing. 
We simulated the HYP1 action on temporal lattices with N t = 4 and N t = 6 time extensions. 
The change in lattice spacing corresponds to about a factor of 5 change in lattice volumes. 
At present we cannot match the physical scales between the N t = 4 and N t — 6 simulations. 
Nevertheless, comparing temperature regions where similar changes in the observables occur 
(e.g. the region where the chiral condensate starts increasing as the temperature is lowered) 
we notice that we can update about 5 times as many links on the N t = Q lattices than on 
the N t = 4 lattices. This indicates that the physical volume of the updated region remains 
constant as the continuum limit is approached. 

III. FINITE TEMPERATURE RESULTS 

In this finite temperature study we performed simulations at Nt — 4 on 8 3 x 4, 10 3 x 4 
and 16 3 x 4 lattices and at Nt = 6 on 16 3 x 6 lattices. At Nt = 4 we considered both the 
APE1 and HYP1 actions and we will compare our results to thin link simulations. We used 
the data from quenched spectroscopy simulations shown in figure |1] to match approximately 
the quark masses of the different actions. The bare mass values used in our simulations are 
listed in table [E|. At the lowest quark mass value we simulated only the HYP1 action as 
it became computationally very demanding to simulate both the thin and APE1 actions. 
At Nt — 6 we simulated only the HYP1 action with a mass that corresponds to the lowest 
mass at N t = 4. To get a feeling for the physical value of the masses, assuming a critical 
temperature T c ~ 150 — 170 MeV for four flavor QCD [§] and using Z = 1 for the mass 
renormalization factor we estimate that the lowest quark mass simulations corresponds to 
m q ~ 6 — 7 MeV. Even though this estimate can easily have a factor of two to four error, 
this value is still very close to the physical light mass values. 

The observables that we measure in our simulations are the average Polyakov loop 



L = pE Tr IT u (n,t),o, 



N t -l 



(10) 



the chiral condensate 



1 



1 



10 



ijjijj 



NfN t 



tr M-\V) 



lONfN t 



Y^R\M-\V)R t 



(11) 



* i=i 
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Action/Set 


Thin 


APE1 


HYP1 


N s 


Set 1 (N t = 4) 


0.06 


0.08 


0.1 


8 


Set 2 {N t = 4) 


0.024 


(0.032) 


0.04 


8, 10 


Set 3 {N t = 4) 


(0.006) 


(0.008) 


0.01 


8, 10, 16 


Set 3 {N t = 6) 


(0.004) 


(0.0053) 


0.0067 


16 



Table II: The bare quark mass values am q for the different actions that correspond to approximately 
the same physical mass. An entry in parenthesis indicates that we did not simulate that action in 
the data set. In the last column we list the spatial extension N s of the lattices in our simulations. 
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Figure 2: Thin link action simulated on 8 3 x 4 lattices at quark mass am q = 0.06 (Set 1 in table 
U) • The real part of the Polyakov line and the chiral condensate < tp ip > are plotted as a function 
of p. 



which we evaluate using 10 Gaussian random vectors Ri per configuration and the discon- 
nected chiral susceptibility 

Xw,disc = ^N*N t {< {^f > - < > 2 ). 
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Figure 3: APE1 action simulated on 8 3 x 4 lattices at quark mass am q = 0.08 (Set 1 in table 0). 
The real part of the Polyakov line and the chiral condensate < iptf) > are plotted as a function of 
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Figure 4: HYP1 action simulated on 8 3 x 4 lattices at quark mass am q = 0.1 (Set 1 in table 0). 
The real part of the Polyakov line and the chiral condensate < ipip > are plotted as a function of 
(5. Dotted symbols indicate data points where the lattice spacing is too large to expect continuum 
like behavior. 
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Action 


P 


am q 


r /a 


y/aa 


a [fm] 


A. 


m G /m p 


HYP1 


5.0 


0.1 




0.51(3) 


0.23(2) 


0.144(6) 


0.66(2) 


HYP1 


5.2 


0.1 


3.03(4) 


0.37(1) 


0.16 


0.057(6) 


0.719(8) 


thin 


5.2 


0.06 


2.85(2) 


0.40(1) 


0.18 


0.63(3) 


0.56(2) 



Table III: Physical scale and flavor symmetry violation of the dynamical HYP1 action. For com- 
parison we list the corresponding values for the thin link action whose bare parameters are chosen 
to approximately match the physical scale and Goldstone pion mass of the (3 = 5.2, am q = 0.1 
HYP1 simulations. 

A. Phase structure and flavor symmetry 

Let us look first at the largest mass set, Set 1 in table [E|. Figures ||, || and [| show the 
real part of the Polyakov line and < tptp > for the three actions on 8 3 x 4 lattices. The thin 
link action shows a pronounced discontinuity at j3 ~ 5.1. The jump in the Polyakov line 
is about 0.4, in < ipip > about 0.2. The APE1 action also gives indication of a first order 
phase transition at /3 ~ 4.85 but the discontinuity of the Polyakov line is only about 0.2, of 
< ipip > about 0.1. Also, the discontinuity is at a lower (3 value with respect to the thin link 
action. Although we do not have a scale determination for the APE1 action, it is unlikely 
that the physical value of the phase transition temperature is unchanged. As figure |] shows 
the discontinuity completely disappears with the HYP1 action. Both the deconfmement 
and chiral order parameters change smoothly from the low to the high temperature phase. 
Moreover the susceptibility X^,disc does not show a peak as a function of the coupling (3 in 
the investigated range. 

The observation that the phase transition weakens as the gauge connections between 
fermions become smoother and eventually disappears for our smoothest HYP1 action is a 
strong indication that the phase transition of the thin link action is a lattice artifact. To 
support this further we determined the physical scale of the HYP1 action at two coupling 
values, (3 = 5.0 and /3 = 5.2, with mass am q = 0.1 on 8 3 x 24 lattices. On the same 
configurations we also measured the meson spectrum. Our results are summarized in table 
PTI] where we list the values for the Sommer scale r /a, the string tension ay/a, the relative 
mass splitting A^ = (m^ — m G )/m G between the lightest non-Goldstone pion 71^5 and the 
Goldstone pion G and the pion to rho mass ratio. At (3 = 5.0 the violations of rotational 
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symmetry in the static potential are substantial and we were not able to determine reliably 
r /a. We could obtain a value for the string tension though from which we get the lattice 
spacing quoted in table |T|. At (3 = 5.2 both quantities r Q /a and ay/a give a consistent 
value for the lattice spacing. For comparison in table [Til] we also give the corresponding 
values for the thin link action at (3 — 5.2, am q = 0.06. This quark mass value corresponds 
approximately to the am q = 0.1 HYP bare mass. As the results in table |TTT| show the scale 
of the HYP and thin link actions are fairly close at the same gauge coupling once the bare 
quark mass values are matched. 

Again, we see that the flavor symmetry restoration predicted by quenched simulations 
agrees with the dynamical results. The quenched study of Ref. |16|] predicted, at a lattice 
spacing a ~ 0.17 fm, a mass splitting = 0.049(8) with the HYP1 valence action at 



am q = 0.1. As table |I| shows the dynamical simulations confirm this result. Flavor 
symmetry with the HYP1 action is an order of magnitude better than with the thin link 
action. 

The lattice spacing of the HYP action at f3 — 5.0, am q = 0.1 is already 0.23fm, 1.4 times 
larger than the lattice spacing at (3 — 5.2, am q = 0.1. Naive scaling implies that the lattice 
spacing is above a = 0.3fm at (3 — 4.8. Considering even lower (3 values does not make 
much sense physically. Once the lattice spacing is comparable or larger than the typical 
instantons on the configuration one cannot expect to see continuum like chiral behavior. We 
used dotted symbols in figure |] to plot the (3 < 4.8 results to indicate that those data points 
are not likely to show continuum behavior. 



B. The phase diagram with the HYP action 

The HYP1 action shows only a crossover with am q = 0.1 quark mass. This might not be 
that surprising if we look at the endpoints of the first order phase transition lines for the 
thin link action quoted in Qj. At N t = 4, am max ~ 0.073 while at N t = 6, am max ~ 0.021, 
showing a substantial scale breaking. The physical mass of Set 1 is lower than m max at 
N t = 4 but larger than m max one would predict from the N t = 6 simulations assuming 
scaling. 

In this section we consider the HYP1 action at smaller quark masses to see if we can 
find a first oder phase transition. The quark mass of Set 2 is, in physical units, below the 
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Figure 5: HYP1 action simulated on 8 3 x 4 lattices at quark mass am q = 0.04 (Set 2 in table ||). 
The real part of the Polyakov line and the chiral condensate < ipip > are plotted as a function of 
(5. Dotted symbols indicate data points where the lattice spacing is too large to expect continuum 
like behavior. 
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Figure 6: HYP1 action simulated on 8 3 x 4 (crosses), 10 3 x 4(diamonds) and 16 3 x 4 (filled octagons) 
lattices at quark mass am q = 0.01 (Set 3 in table |ll|). The real part of the Polyakov line and the 
chiral condensate < ipip > are plotted as a function of /?. Dotted symbols indicate data points 
where the lattice spacing is too large to expect continuum like behavior. 
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observed endpoint m max of the N t = 6 thin link transition. Numerical results for the thin 
link action at am q = 0.024 on 8 3 x 4 lattices show a strong first order phase transition. Yet 
the results with the HYP1 action on 8 3 x 4 lattices plotted in figure [5| show only a very 
broad crossover. 

Looking for the phase transition signal we simulated the HYP1 action at the even lower 
quark mass of Set 3, am q = 0.01. At this low value of the quark mass finite volume effects 
could be important. For example the first order phase transition observed with the thin link 
action at N t = 6 on 16 3 x 6 lattices at a quark mass am q = 0.01 in Ref. || washes away 
if we repeat the simulations on 8 3 x 6 lattices. To make sure that our results are not due 
to finite size effects we simulated the HYP1 action at quark mass am q = 0.01 on 8 3 x 4, 
10 3 x 4 and 16 3 x 4 lattices. The results are shown in figure ||. There is no sign of significant 
deviations as the spatial volume is increased by a factor of eight and the phase diagram 
shows a very broad crossover. We conclude that there is no signal for a first order phase 
transition at N t = 4 for quark masses down to approximately the physical light quark mass 
in the physically relevant gauge coupling range. 

Since we found no sign of a phase transition, only a gradual change in the Polyakov loop 
and chiral condensate expectation values, one might wonder if we found a chirally broken 
phase at all. To study this question we consider the chiral condensate at fixed gauge coupling 
as the function of the bare quark mass. For small quark mass the condensate is expected 
to be linear in the mass suggesting to look at the quantity < t/jip > /m q as the function of 
l/m q 

S±t2L = lL + c (12) 
m q m q 

where £ is the value of the condensate extrapolated to zero quark mass. Figure [7] shows the 
possible situations one might encounter. The solid line in the figure separates the chirally 
broken (low temperature) and chirally symmetric (high temperature) phases. This line 
could describe a phase transition or only a crossover. The points at temperature T 3 are 
in the chirally symmetric phase, S =< ^pip > \ mq =o = 0, the linear form of eq.(|i~2"D should 
predict a constant value for < ipip > l m q- The points at 7\ are all in the chirally broken 
phase, £ is finite and positive, eq.(|i~2|) should predict a linear dependence for < ipip > /m q 
with a positive slope. The case at T 2 is more complicated. The points below the solid line 
(open symbols) are in the chirally symmetric phase and should behave the same as the open 
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Figure 7: A schematic phase diagram in the temperature - quark mass plane. Open symbols are 
in the chirally symmetric, filled symbols are in the chirally broken phase. The solid line represents 
the phase transition or crossover between the phases. 

symbol points at T 3 giving £ = in the extrapolation. The points above the solid line (filled 
symbols) are in the chirally broken phase and the linear form eq.(|l2D should describe the 
mass dependence in the chirally broken phase with a positive condensate at the transition 
point. The quantity £ in eq.(|l2|) however is not this quantity but its linearly extrapolated 
value to zero quark mass. £ could therefore take any value, even negative ones, if at the 
transition point the condensate is small or the quark mass is large. 

We have analyzed the chiral condensate on the 8 3 x 4 configurations according to eq.(p^). 
In addition to the quark mass values we presented above we simulated the system at a few 
additional ones to make the extrapolation meaningful. Figure |8] shows < Tpip > /am q as 
the function of l/am q at (3 = 4.8 , /3 = 4.9 and (5 = 5.0. These (3 values correspond to the 
situations described at T 2 , T 2 and T\, respectively, in figure |7|. For easier identification we 
have used the same symbols in the two figures. At (3 — 4.8 and 4.9 the slope £ is negative 
for am q > 0.025(0.03), and £ is consistent with zero for am q < 0.025(0.03) indicating that 
there is a phase boundary separating the chirally broken and symmetric phases around 
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Figure 8: < ipip > /m q as the function of l/m g at f3 = 4.8, /3 = 4.9 and /? = 5.0. All data are on 
8 3 4 lattices. 



am n 



0.025(0.03). We have found negative S slope at all /3 < 4.9 coupling though the 



errors at (3 < 4.7 are too large to distinguish E from zero. At (3 = 5.0 all data points 
are consistent with a vanishing slope, S = 0, suggesting that all points are in the chirally 
symmetric phase. We have found similar behavior for (3 > 5.0. These observations suggest 
that the simulations have found the chirally broken phase with all investigated quark mass 
values. They also suggest that the chiral condensate at the crossover is small and that the 
chiral am q = phase transition, if exists, is at a fairly low (3 value on N t = 4 lattices. These 



results are consistent with the interacting instanton liquid model predictions ||TTJ, |12| but 
not with the thin link action lattice results. 

In order to understand the phase structure of four flavor fermions it would be important 
to determine the temperature at the am q = phase transition point and locate the phase 
transition line at finite quark mass. This, however, should not be done on N t = 4 lattices. 
As we pointed out in the introduction, simulations at Nt — 4 might be questionable even 
with chirally symmetric actions because the lattice spacing is too large to support the correct 
topological structure of the vacuum. Furthermore, decreasing the quark mass while staying 
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Figure 9: HYP1 action simulated on 16 3 x 6 lattices at quark mass am q = 0.0067 (Set 3 in table 
U) . The real part of the Polyakov line and the chiral condensate < Tptp > are plotted as a function 
of (3. 

at N t = 4 makes flavor symmetry violations larger even with the HYP1 action. We believe 
that the disappearance of the first order phase transition is due to the improved flavor 
symmetry of the HYP1 action. These considerations indicate that we have to consider 
lattices with N t > 6 temporal extension. 

In figure || we show our first results for N t = 6 thermodynamics with HYP1 action on 
16 3 x 6 lattices at a quark mass am q = 0.0067. This mass matches approximately the lightest 
quark mass of Set 3 at N t = 4. The critical region, which we identify as the region where 
< ipip > starts increasing and where the autocorrelation times are the largest, moves to 
higher (3 values than at N t = 4. Again, there is no sign of a first order phase transition, the 
Polyakov line and the chiral condensate are smooth, indicating a crossover as observed at 
N t = 4. The susceptibility Xiptp,disc does not show any peak in the investigated range. With 
a shift A/? rs 0.2 the susceptibility agrees with the one obtained on the N t = 4 lattices. Also 
the chiral condensate normalized by the quark mass, < Tpip > /m q matches the N t = 4, 



0.01 data if the gauge coupling is shifted by A/3 = 0.2. This shift in the gauge 



coupling for a scale change of | = 1.5 is consistent with what we observed between the 
(3 = 5.0 and j3 = 5.2, am q = 0.1 zero temperature simulations. 

Our N t — 6 results are only preliminary. Obviously the simulations should be repeated at 
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different quark mass values in order to do the extrapolation according to eq. (|I"2"D but that 
is beyond the scope of this paper. Our goal here was to demonstrate the striking difference 
between the thin link and HYP actions. 

IV. SUMMARY 

In this paper we studied the finite temperature phase diagram of four flavor fat link 
staggered fermion actions on lattices with N t = 4 and N t = 6 temporal extents. We used a 
new smearing, the hypercubic blocking (HYP), for the fat links. Hypercubic blocking mixes 
the gauge links only within hypercubes attached to the original link minimizing thereby the 
lattice artifacts arising from extended smearing. The action with one level of hypercubic 
blocking (HYP1) considerably reduces the flavor symmetry breaking of the standard thin 
link action but does not change the local properties of the gauge fields beyond distances 
r/a > 2. We presented a new algorithm for simulating the HYP1 action. This algorithm is 
very simple because it is based on the overrelaxation and heatbath algorithms for the pure 
gauge action. It is efficient as the quark mass is lowered allowing simulations which would 
be impractical with thin link action. 

We found that improved flavor symmetry, presumably the presence of 15 near-degenerate 
Goldstone bosons, changes the finite temperature phase diagram of four flavor QCD. Our 
simulations with the HYP1 action do not show any sign of a first order phase transition 
even at physical quark masses of a few MeV. We found indications for a broad crossover- 
like transition between the chirally broken and symmetric phases. The phase transition 
temperature at zero quark mass is too small and cannot be identified on N t = 4 lattices 
even with chirally improved actions. 

Our algorithm to simulate the HYP1 action can be generalized to arbitrary number of 
flavors by using a polynomial approximation for the fermionic action ||19|| . 
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